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Q , The early X-ray afterglow for a significant number of gamma-ray bursts detected 

^ ! by the Swift satellite is observed to have a phase of very slow flux decline with time 

^ : (F^oc r'^ with 0.2 < a < 0.8) for 10^^ s < ? < 10^ s, while the subsequent decline 

^ ! is the usual 1 < as < 1.5 behavior, that was seen in the pre-Swift era. We show that 

^ i this behavior is a natural consequence of a small spread in the Lorentz factor of the 

[ ejecta, by a factor of ~ 2 - 4, where the slower ejecta gradually catch-up with the 

5 I shocked external medium, thus increasing the energy of forward shock and delaying 

^ I its deceleration. The end of the "shallow" flux decay stage marks the beginning of the 

^ I Blandford-McKee self similar external shock evolution. This suggests that most of the 

^ I energy in the relativistic outflow is in material with a Lorentz factor of ~ 30-50. 

O I Subject headings: gamma-rays: bursts — shock waves — hydrodynamics 

+^ ■ 

> ' 

. 1. Introduction 

_C3 I Among the discoveries made by the Swift satellite within a few months of its launch is the 

observation that a fraction of long duration gamma-ray bursts (GRBs) go through an early phase 
of relatively slow decline in the X-ray afterglow flux that typically starts at a few minutes after the 
burst and lasts for about an hour (Nousek et al. 2005). This phase is followed by a somewhat faster 
and more typical flux decay that satisfies the expected relation between the temporal decline index 
a and the spectral index P, where oc z/"^r", similar to what was observed before the Swift era 
when the monitoring of the afterglow light curves started at least several hours after the ORB. The 
spectral index does not seem to undergo any change when the light-curve transitions (at ?break,2 ~ 
10"^ s) from a shallow decline (^2) to the "regular" decline (03). It has been argued convincingly 
by a number of authors that the more slowly declining lightcurve, like the "regular" flux decay rate 
that follows it, are both produced by the shock heated circum-burst medium (Nousek et al. 2005; 
Panaitescu et al. 2005; Zhang et al. 2005). The shallow X-ray flux decay is widely attributed to 
energy injection into the afterglow shock, which may be caused by either a long lived activity of 
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the central source, or a short lived central explosion that produces ejecta with some distribution of 
Lorentz factor (LF), cf. (Nousek et al. 2005; Panaitescu et al. 2005; Zhang et al. 2005). In either 
of these scenarios the deceleration of the afterglow shock is reduced due to the energy being added 
to it, and this in turn produces a slowly declining light curve. 

Slow decline of the early optical lightcurve was also reported before the Swift era, although 
quite rarely, e.g. GRB 021004 (Fox et al. 2003). Li & Chevalier (2003) have made a good case that 
the slowly declining early optical lightcurve of GRB 021004 could naturally occur for a stellar wind 
type external medium when the optical band is below the characteristic synchrotron frequency (and 
above the cooling break frequency). However, this explanation is unlikely to work for the slowly 
declining X-ray lightcurves since the X-ray band at ~ 1 hr after a GRB is expected to lie well 
above the synchrotron characteristic frequency, and furthermore there is no evidence for a change 
in the spectral slope across the break in the X-ray lightcurve at ?break,2- Fox et al. (2003) argued that 
the early flat optical lightcurve of GRB 021004 is either due to energy injection into the afterglow 
shock, or due to angular inhomogeneity ("patchy shell"; Kumar & Piran 2000). The latter was 
favored as it also provided a good explanation for the fluctuations (or "bumps") that appear later 
in the optical lightcurve of GRB 021004 (see also Nakar, Piran & Granot 2003). However, the 
sparse early afterglow data in the pre-Swift era made it difficult to distinguish between the different 
explanations. 

A long lived activity of the central source is not very appealing since it would require the 
source to be active up to several hours after the GRB, with a very smooth temporal behavior, where 
most of the energy is in the outflow that is ejected around ?break.2 ~ lO'* s (see however Dai 2004); 
this makes the problem of the observed high efficiency for converting kinetic energy to gamma-ray 
radiation much worse (Nousek et al. 2005). Another interesting way to produce an early flat phase 
in the afterglow light curve (Eichler & Granot 2005) is by a line of sight that is slightly outside the 
(sharp) edge of a roughly uniform jet (Granot et al. 2002; Granot, Ramirez-Ruiz & Pema 2005). 
This would, however, naturally be accompanied by a weaker and softer prompt emission, perhaps 
resulting in an X-ray flash or X-ray rich GRB rather than a classical GRB; the more pronounced 
this effect is the flatter and longer lived the slow X-ray afterglow decay phase should be. Initial 
inspection of the data does not show such a correlation, suggesting that viewing angle effects are 
probably not the predominant cause of the early slow decay phase in the X-ray afterglows, at least 
under the simplest assumptions.^ 

It is natural to expect that matter ejected in any explosion will have a range of velocities or 



^Eichler & Granot (2005) point out that viewing angle effects might still be the dominant cause of the flat early 
decay of the afterglow light curves if along some lines of sight the kinetic energy in the afterglow shock is very low 
while the energy in gamma-rays remains high. 
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LFs. After a while (on a time scale, in the observer frame, of order a few times the duration of the 
central engine activity) the ejecta will rearrange themselves such that the fastest moving plasma is 
at the head of the outflow and the slowest at the tail end. This can occur either through internal 
shocks within the outflow, or by a smooth decrease in the LF of the outflow toward the end of the 
central source activity. If the ejecta have a finite range of LPs, the slower ejecta would gradually 
catch up with the shocked external medium, injecting energy into the forward shock. If the slower 
ejecta carry more energy than the faster ejecta, then this added energy would gradually increase 
the energy of the afterglow shock, causing it to decelerate more gradually. Once the energy in 
the lower LF ejecta becomes small compared to the energy already in the afterglow shock, the 
blast wave evolution becomes impulsive (i.e. the subsequent small amount of energy injection 
hardly effects the evolution of the forward shock), and if radiative losses are unimportant then it 
approaches the adiabatic Blandford & McKee (1976) self-similar solution. This occurs when the 
LF of the afterglow shock drops slightly below Fpeak, the LF where dE/d\nT peaks and where 
most of the energy in the outflow resides. 

In this paper we use the Swift data to determine the time dependence of the blast wave LF. 
We find that the LF typically drops by a factor of ~ 2-4 during the shallow decline phase. This 
is consistent with the basic picture suggested above, where a finite LF distribution for the ejecta 
causes a more gradual decline of the forward-shock LF, which gives rise to a shallow light-curve, 
and is an intermediate transition stage before the onset of the adiabatic Blandford-McKee solution. 

2. Dependence of Burst Kinetic Energy on Lorentz Factor 

The emission from an external shock can be described in terms of the shock front LF (F) 
and the density profile of the circum- stellar medium (CSM). For a uniform CSM the synchrotron 
characteristic frequency {u,„), the cooling frequency {Uc) and the flux at the peak of the spectrum 
(^i/,max), in the observer frame, are proportional to F'*, T~^t~^ and ^^F*^ respectively, where t is 
the observed time. The flux at a frequency between the Vm and Vc is proportional to t^T^'^P and 
for the observed band above v„ and z/^ the flux scales as t^T^'^P. The observed flux is strongly 
dependent on F and therefore even a small deviation from the F oc t~^l^ scaling has a very large 
effect on the observed light-curve. The observed flux has a weaker dependence on F for a wind 
like density stratification of the CSM; the flux in the two regimes considered above scales roughly 
as T^+Pt^^-p)n and v^+Pt^P-^V^, respectively. 

More generally, for a power law external density profile, pext = Ar""*, we have /v,max oc TBE?~^ oc 
p2^3-3V2 ^ Y^-3kp-3k/2^ ^ p^^2 ^ ^^j^-k/i ^ Y^^-^l^, 7, OC 1 /VBH and Vc oc TB^l oc F-^fi-^r^ oc 
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R^k/i-p-Af-i ^ p3*-4^-2+3V2 Therefore, 



(1) 



Assuming that the LF distribution for the ejecta is E{> V) oc V~'^, we find g = -dlogV / dlogt 
is smaller by an amount 6 compared to the standard value of (3-k)/2(4-k), i.e. 3/8 (1/4) for a 
uniform (wind) CSM, where 



S = 



0-k)a r3«/[8(8+a)] ^ = 0, 

2(4-.)[2(4-.).a]-|^/^^^^^^^^ ..2. 



(2) 



The deviation to the LC temporal power-law index (Aa) from the standard case of Blandford- 
McKee self-similar solution {a) is easily related to 5. For v,n<^ <^c^^ have 



Aq;= 6- 



while for v > max(i^;„, u^). 



Aa^ 4-k+ 



p(4-k)-5k 



p(4-k) 



S = 



■3(3+/7)a/[4(8+fl)] k = 0, 
(l+p)a/[4(4+fl)] k = 2, 



•3(2+/7)a/[4(8+fl)] k = 0, 
(2+p)a/[4(4 + a)] k^2. 



(3) 



(4) 



We next calculate S for a number of Swift detected GRBs with a shallow LC using the ob- 
served spectral index and the change in the temporal power-law index for the X-ray lightcurve 
(Aa) between the shallow and the "regular" parts of the LC. The results for S, and the change to 
the LF during the shallow LC are shown in Table 1. 

We note that the change to F during the shallow phase of the LC was calculated using the 
appropriate dependence of F on for a uniform CSM this is t^^^^^. It can be seen in Table 1 that F 
changes by a factor ~ 2-4 for all the bursts, with a uniform CSM, during the shallow LC phase; 
these numbers change only by a small amount even if we take the forward shock emission, and the 
shallow decline, to begin at the end of the GRB. 

The function dE/dlnT peaks at Fpeak ~ r(?break,2), the LF of the forward shock at the end of 
the shallow decline phase of the X-ray lightcurve. For F > Fpeak, dE/dlnV oc T~". The power-law 
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index a is given in table 1 for a number of bursts detected by Swift and lies between ~ 1 and ~ 2.5 

if the CSM has uniform density (Nousek et al. 2005, report similar values - 5- 1 in their notation); 
a>5 if the medium in the vicinity of GRB is taken to be a wind-CSM'' or alternatively the central 
source has to be active for several hours with little variability and a roughly constant rate of energy 
output in relativistic outflow - neither of these possibilities seem very plausible and so the case of 
a wind-CSM is not considered any further in this paper. For F < Fpeak the function dE/dlnT (xT'^ 
should decrease with decreasing F (i.e. ft > 0) as otherwise slower moving ejecta will continue 
to add substantial amount of energy to the forward shock thereby retarding its deceleration and 
slowing down the decline of the lightcurve. Since the spectral index and the lightcurve power- 
law decay index after the end of the shallow decline phase obey the relationship expected for an 
adiabatic forward shock evolution we conclude that indeed b> 0, but its exact value is otherwise 
unconstrained. Radio calorimetry for a number of GRBs has concluded that there is not a whole 
lot of energy in GRBs in the form of mildly relativistic ejecta with F ~ 2 (e.g. Berger et al. 2004; 
Frail et al. 2005). This further strengthens our conclusion that b > 0, and that this scaling might 
extends to F ^ 2. We note that for a given total energy in the explosion of order 10^^ erg the 
relation dE/dlnV oc F"", with a ~ 1.5, must turnover at some F of order 10 or so otherwise the 
energy in the relativistic ejecta will exceed the total available energy (energy in relativistic ejecta 
with F > Fpeak ~ 50 is of order 10^^ erg). We have now considerable body of evidence that long 
duration GRBs are accompanied by a supernova of Type Ic, which expels a few solar masses of 
material at velocities of order 10"^ km s"'. Thus, dE/dlnu, where u = (3T = (F^- 1)^/^, must again 
turnover over and have a peak at m ~ 0.05. Putting all these together we show a schematic behavior 
of E(l3T) (i.e. dE/d\nu) in Figure 1. 



3. Conclusion 

We have pieced together the distribution of energy in gamma-ray burst ejecta as a function of 
the four- velocity a = /5F = (F^- 1)^/^ for 0.1 < m < 10^. The distribution function, dE/dlnu, has 
two peaks: one at m ~ 0.1 and another at m ~ 30-50. For u > 50, it falls off as dE/dlnu oc m"'' 
with a ~ 1 - 2, as is determined from the shallow decline of the X-ray lightcurve at early times 
(10^^ s < ? < 10"* s) observed for a good fraction of bursts detected by the Swift satellite. The 
distribution at low w ~ 0. 1 is obtained by observations of supernovae Ic that are associated with 
GRBs. In the intermediate regime of 1 < m < 30 the shape of the distribution function dE/dlnu is 
very uncertain, but we argue that it is likely to be at least flat or slowly rising in this range. 



''The total amount of energy injection during the shallow dechne phase is independent of the stratification of the 
circum-stellar medium. However, for a wind-hke CSM a>5 whichmeans that t/E/t/lnr must have a very narrow peak 
of width (5hir <C 1 which is unhkely to be reahzed in nature. 



-6- 



A prediction of this model for the shallow decline of X-ray light-curve, that is based on a 
deviation (at early times) from the constant energy Blandford-McKee self-similar solution, is that 
we should see a roughly similar shallow decline in the optical band over the same time interval as 
in the X-ray data. Since the optical and the X-ray bands in general lie in different segments of the 
synchrotron spectrum, and because the energy added to the forward shock by slower moving ejecta 
should be accompanied by a mildly relativistic reverse shock that could provide some added flux to 
the optical lightcurve, the rate of decline for the X-ray and optical lightcurves should be similar but 
not identical in this model.'^ A highly magnetized outflow could significantly weaken the reverse 
shock (or even eliminate it altogether) and thus suppress its emission. The alternative explanation 
of a viewing angle slightly outside the edge of the jet would lead to a gradual steepening of the 
afterglow lightcurve (i.e. a gradual increase in a, as the beaming cone of the afterglow emission 
gradually approaches and eventually encompasses the line of sight), while a steeper break in the 
light curve is possible (and arguably, might also be expected) in the model described in this work, 
when the stage of energy injection into the afterglow shock ends. 

The challenge posed for GRB/SNe models is to understand what physical processes give 
rise to a ~ 2 and why the LF distribution of the ejecta peaks at a value roughly Fpeak ~ 30 - 50. 
Understanding these results should help illuminate the processes operating during the period in 
which the central engine of gamma-ray burst is active and the interaction of the relativistic outflow 
with the collapsing star and its immediate surroundings. 
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Table 1 . The Change in the Lorentz Factor During the Energy Injection Episode 



GRB# 


790/s 


Wak,l/s 


*break,2/s 


a 


5 


^ min ^ 


^max 


-Tpeak * 


To* 


050128 


13.8 


<230 


17205^» 


1.1 ±0.2 


0.045 ±0.07 


2.0 


4.9 






050315 


96.0 


400 ±20 


12000 ±400 


2.4 ± 0.1 


0.094 ±0.03 


2.6 


3.9 


30 


117 


050319 


10.0 


370 ±15 


40000 ±300 


1.6 ±0.5 


0.063 ±0.016 


4.3 


13.3 


21 


279 


050401 


33.0 


< 127 


5500+1150 

-'■'""-1050 


1.7±0.1 


0.066 ± 0.003 


3.2 


4.9 


58 


284 


050416a 


2.4 


< 80 


1 350^20™ 


1 . 1 ± 0. 1 


0.043 ±0.004 


2.6 


8.2 


28 


230 


050607 


26.5 


510 = 50 


6400 ±400 


1.5 ±0.1 


0.059 = 0.004 


1 ~> 


5.7 







Note. — The relevant data were taken from Nousek et al. (2005). AU of the calculated quantities reported in this table - 
a, 5, ^ and Fq - assume a uniform density medium in the vicinity of these bursts, t Here ^ is the ratio of the Lorentz factor 
of the afterglow shock at the start and at the end of the "shallow part" of the X-ray light-curve, and its value is estimated to 
be between ^^n = (Wak,2Ahreak,i)^^^^"*"' and fmax = (Wak,2/79o)'^***°'; * These values of Tpeak = r(fbreak,2) were estimated 
only for the GRBs with known redshifts, by using equation 9 of Nousek et al. (2005) where the isotropic equivalent kinetic 
energy at ?break.2 was taken to be equal to iJ-y^so. and the external density was taken to be n = 1 cm~^; * The initial Lorentz 
factor is simply estimated by Fq = SmaxFpeak- 
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Fig. 1. — Schematic figure showing dE/d\nu, in units of 10 erg, as a function of m = F/? = 
(F^- 1)^/^. It has one relativistic component (soUd line) with total energy ~ 10^* erg and peak at 
M ~ 30-50 that produces the gamma-ray burst and the afterglow radiations. The power-law index 
above the peak for this component is well constrained by the X-ray data (the shallow part of the 
light-curve) and is ~ -1.5 (see table 1). The slope below the peak is not constrained and is taken 
to be 1 ; in reality it can be close to zero, as the only constraint we have is from late time radio 
afterglow observations which suggests that there is not a lot of extra energy in material moving 
with Lorentz factor of order 2. The second component (dashed curve) shows schematically the 
kinetic energy in non-relativistic ejecta in the supernova accompanying the GRB; the peak for this 
component is taken to be ~ 10"^ km s~\ the typical velocity for SNe Ic ejecta, and the energy is 
~ 10^2 erg. 



